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Abstract 

We consider an electron moving in a periodic potential and subject to 
an additional slowly varying external electrostatic potential, ^(ea::), and 
vector potential A{ex), with a; £ R"* and e <^ 1. We prove that associ- 
ated to an isolated family of Bloch bands there exists an almost invariant 
subspace of L^CMJ^) and an effective Hamiltonian governing the evolution 
inside this subspace to all orders in e. To leading order the effective Hamil- 
tonian is given through the Peierls substitution. We explicitly compute 
the first order correction. From a semiclassical analysis of this effective 
quantum Hamiltonian we establish the first order correction to the stan- 
dard semiclassical model of solid state physics. 
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1 Introduction 



A central problem of solid state physics is to understand the motion of electrons 
in the periodic potential generated by the ionic cores. While the dynamics is 
quantum mechanical, many electronic properties of solids can be understood 
already in the semiclassical approximation |AsMel IKoL IZa| . One argues that for 
suitable wave packets, which are spread over many lattice spacings, the main 
effect of a periodic potential Vr on the electron dynamics consists in changing the 
dispersion relation from the free kinetic energy Eh-cdk) = ^ fc^ to the modified 
kinetic energy En{k) given by the n^^ Bloch band. Otherwise the electron 
responds to slowly varying external potentials A, (/) as in the case of a vanishing 
periodic potential. Therefore the semiclassical equations of motion read 

r = \/En{K), k = -V(t){r) +r X B{r) , (1) 

where r G R'^ is the position of the electron, k — k — A{r) its kinetic momentum 
with k its Bloch momentum, —Vcj) the external electric field, and B — V x A 
the external magnetic field. Note that there is a semiclassical evolution for each 
Bloch band separately. (We choose units in which the Planck constant h, the 
speed c of light, and the mass m of the electron are equal to one, and absorb 
the charge e into the potentials). 

One goal of this article is to understand on a mathematical level how these 
semiclassical equations emerge from the underlying Schrodinger equation 

iedt^p{x,t) ^ (i(-iV^ - A(ex))Vyr(a;) +(/)(ex)) 

= H'^{x,t) (2) 

in the limit e ^ at leading order. Here the potential Vr : ^ K is peri- 
odic with respect to some regular lattice F. F is generated through the basis 
{71,72,73}, 7j G M^, i.e. 

r = |a;e]R'^:a; = X]j=i'^j 7j ^'^^ some a G Z'^j , 

and Vr(a; + 7) = Vr{x) for all 7 G F, x G R'^. The spacing of the lattice F defines 
the microscopic spatial scale. The external potentials A{ex) and (j){ex), with 
A : ^ R'^ and </) : R'^ ^ R, are slowly varying on the scale of the lattice, as 
expressed through the dimensionless scale parameter e, e ^ 1. In particular, 
this means that the external fields are weak compared to the fields generated 
by the ionic cores, a condition which is satisfied for real metals even for the 
strongest external electrostatic fields available and for a wide range of magnetic 
fields, see |AsMe| . Chapter 12. 

In solid state physics the derivation of the semiclassical model ^ received a 
lot of attention during the 1950s to the 1970s. We mention representatively the 
work by Luttinger |Lu| . Kohn |Ko| . Blount Bli", 'BI2' and Zak |Za| . As late as 
1962 Wannier |Wa| argues that the derivation of from jSJ is still incomplete. 

On the mathematical side the semiclassical asymptotics of the spectrum 
of have been studied in great detail by Gerard, Martinez and Sjostrand 
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|GMS| with predecessors 'BcRa' 'Bui |HeSj[ INe| . The large time asymptotics 
of the solutions to ^ without external potentials is studied in jAsKn| and 
the scattering theory is developed in GeNi,. However for the dynamics of 
wave functions, our interest here, the results are modest. In |GMMP| the case 
(j) = 0, A = is considered, in 'HST' and |BMP| a proof is given for A = 0, 
which leaves out many interesting applications. The method of Gaussian beams 
is developed in ('■HI" for a weak uniform magnetic field and in |DGR| for 
magnetic Bloch bands. 

In fact, as our title indicates, we are more ambitious and plan to derive 
also the first order correction to ^ . The electron acquires then a fc-dependent 
electric moment Anik) and magnetic moment A^„(fc). If the n*^ band is nonde- 
generate (hence isolated) with Bloch eigenfunctions tpn{k,x), the electric dipole 
moment is given by the Berry connection 

A(fc) =i(V'n(fc),VfcVn(fc)) (3) 

and the magnetic moment by the Rammal- Wilkinson term 

M{k)n = \ (VkMk), x(i/p„(fc) - E{k))VkMk)) ■ (4) 

Here (• , •) denotes the inner product in L'^iM.^/T) and iJpcr(fc) is of Q with 
(j) ^ ^ A for fixed Bloch momentum k, see Eq. H17|) . As will be explained in 
detail, the corrected semiclassical equations read 

r = V^f^En^K) - e B{r) ■ MniK)j - s k X ^nin) , 

(5) 

k = -Vr{4>ir) - e B{r) ■ Mn{K)j +r X B{r) 

with riri(fc) = V X An{k) the curvature of the Berry connection. 

The issue of first order corrections to the semiclassical equations of motion 
has been investigated recently by Sundaram and Niu |SuNi| in the context of 
magnetic Bloch bands, see also Chang and Niu |ChNi| . One adds in ||2J) a strong 
uniform magnetic field Bq, i.e. the vector potential Ao{x) — ^Bq x x. If its 
magnetic flux per unit cell is rational, then the Hamiltonian in (0) is still periodic 
at the expense of a larger unit cell and replacing the usual translations by the 
magnetic translations. Equation © remains formally unaltered, only En now 
refers to the energy of the magnetic subband. Instructive plots of ri„ and AAn 
are provided in (SuNi for the particular case of the 2-dimensional Hofstadter 
model at rational flux 1/3. The first order corrections obtained in jSuNi| agree 
with our equation |(SJ), except for the term of order e in the second equation. On 
a technical level magnetic Bloch bands require some extra considerations and 
we defer them to a forthcoming paper IPST3I . 

It has been recognized repeatedly, as e.g. emphasized in |ABL| . that the 
geometric phases appearing in the first order correction contain novel physics as 
compared to the leading order. Bloch electrons are no exception. For example 
for the case of magnetic Bloch bands, the equations of motion Q provide a 
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simple semiclassical explanation of the quantum Hall effect. Let us specialize 
((SJ to two dimensions and take B(r) = 0, (/){r) = —£ ■ r, i.e. a weak driving 
electric field and a strong uniform magnetic field with rational flux. Then, since 
K = k, the equations of motion become r = VkEn{k) + £'^V,n{k), k = £, where 
fin is now scalar, and f-*- is £ rotated by tt/2. We assume initially fc(0) = k 
and a completely filled band, which means to integrate with respect to k over 
the first Brillouin zone M*. Then the average current for band n is given by 

jn=/ dfcf(fc)-/ dk{VkEnik)-£^nnik)) = ~£^ [ dkn^ik). 

Jj^, dk flnik) is the Chern number of the magnetic Bloch bundle and as such an 
integer. Further applications related to the semiclassical first order corrections 
are the anomalous Hall effect JNM] and the thermodynamics of the Hofstadter 
model |GaAv| . 

Our derivation of lO from ||2Jl proceeds in two conceptually and mathemat- 
ically distinct steps. The first step is to obtain an effective Hamiltonian whose 
unitary group closely approximates the solution to the Schrodinger equation ^ 
for e small, in case the initial wave function lies in a subspace corresponding 
to a prescribed family of Bloch bands. Inside the family, band crossings and 
almost crossings are allowed. It is crucial however that for every k the family of 
bands is separated by a gap from the remaining energy bands. Then, associated 
to the given family of bands, there is a subspace H^L^(M^) which is adiabati- 
cally decoupled from its orthogonal complement to all orders in e. The effective 
Hamiltonian generates the approximate time evolution in H^L^(M^). 

Compared to the space-adiabatic perturbation theory developed in |PSTi| , 
as a new element we have to face the fact that the classical phase space is 
(RVr*) X R^, T* the lattice dual to T and R^/T* = M* the first Brillouin 
zone. To come close to the scheme in |PSTi| a natural approach is to use 
the extended zone scheme. Going from one cell to the next, one picks up a 
phase factor which necessitates to generalize the pseudodifferential calculus to 
T-equivariant symbols, see Appendix A. 

The effective Hamiltonian is expanded in an e-independent reference Hilbert 
space. For example, for a nondegenerate band the reference space is L?{M* ,dk) 
and the leading order effective Hamiltonian is given through the Peierls substi- 
tution 

/io(fc,ieVfe) = En{k ~ A(ieVk)) + ^(ieVfc) , (6) 

where iVfc is understood with periodic boundary conditions on AI* . 

The natural second step consists in a semiclassical analysis of the effective 
Hamiltonian. It is a standard result that the unitary group generated by ho is 
well approximated by the semiclassical equations At next order, hQ{k, ieVfc) 
is corrected to ho{k,ieyk) -|- e/ii(A:, ieV^), with hi given in (|^ . However ^ 
is not the semiclassical evolution corresponding to that Hamiltonian. The rea- 
son is that the subspace H'^L^(M^) is mapped to the reference Hilbert space 
L^(M*,dfc) through a unitary operator which itself depends on e. Therefore, 
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the transformation of observables generates an e-dependence in addition to the 
transformation of time-evolved states. If done properly, one arrives at (|5|). 

To give a brief outline of the paper. In Section 2 we discuss the periodic 
Hamiltonian. In particular we recall the unitary Zak transform and state our 
assumptions onVr, A, (j) and the gap condition. In Section 3 we apply the space- 
adiabatic perturbation theory to the present case, using the pseudodifferential 
calculus developed in Appendix A. The semiclassical analysis of the effective 
Hamiltonian including first order is carried out in Section 4. The precise link 
between Q and Q is stated in Theorem^] In Appendix B we show that the 
equations (O are of Hamiltonian form with respect to an appropriate symplectic 
structure. 

2 The periodic Hamiltonian 

In order to formulate our setup we first need to recall several well known facts 
about the periodic Hamiltonian 

i?pcr ■.= -^A + Vr, 

acting in L'^{M.'^), keeping from now on the dimension d arbitrary. The potential 
Vr is periodic with respect to the lattice F. Its dual lattice T* is defined as 
the lattice generated by the dual basis {7J' , . . . , 7^ } determined through the 
conditions 7; ■ 7* = 27r(5y , i,j G {1, . . . ,d}. The centered fundamental domain 
of r is denoted by 

M ^ ^^x eW' : X = E^=i"j 7j for G [-i, i]| , 

and analogously the centered fundamental domain of F* is denoted by M* . In 
solid state physics the set M* is called the first Brillouin zone. In the following 
M* is always equipped with the normalized Lebesgue measure denoted by dk. 
We introduce the notation x — [x] + j for the a.e. unique decomposition of 
a; G K"^ as a sum of [x] E M and 7 G F. We use the same brackets for the 
analogous splitting k = [fc] + 7*. 

We employ a variant of the Bloch-Floquet transform, called the Zak trans- 
form (also Lifshitz-Gelfand-Zak transform). The Zak transform of a function 
^Z- G S{R'') is defined as 

{Uij){k, x) := e-'^"+''^-'=V(a; + 7), {k, x) G (7) 

and one directly reads off from Q the following periodicity properties 

{Ui^){k,y + -f) = {U^,){k,y) for all 7GF, (8) 
(Utp) (fc + 7* , y) = e-'^-'^* {U^) (fc, y) for all 7* G F* . (9) 
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From (|SJ) it follows that, for any fixed k G Mf^, (Utp) {k, •) is a F-periodic function 
and can thus be regarded as an element of L^(T'^), T*^ being the flat torus M'^/F. 
Equation © involves a unitary representation of the group of lattice translations 
on F* (denoted again as F* with a little abuse of notation), given by 

T:F*^iY(i2(T'^)), 7*.^r(7*), (r(7*H(2;) = e'^-^V(2/). 

It will turn out convenient to introduce the Hilbert space 

:={^eiL(K',i'(T')): V'lfc - 7*) = ^(7*) ^(fc)} , (10) 
equipped with the inner product 

(p)h^ = dk (ipik), (y5(fc))i2(T) . 
Jm* 

Notice that if one considers the trivial representation, i.e. r = 1, then Tir is 
simply a space of F*-periodic vector-valued functions over R''. 

Obviously, there is a natural isomorphism between Hr and i^(A'f*, L^(T'')) 
given by restriction from M."^ to M* , and with inverse given by r-equivariant 
continuation, as suggested by @. The reason for working with instead of 
L^(M*, L^(T'')) is twofold. First of aU it allows to apply the pseudodifferential 
calculus as developed in Appendix A. On the other hand it makes statements 
about domains of operators more transparent as we shall see. 

The map defined by {T)) extends to a unitary operator 

U : L2(R'^) -^nr = L^(M\L^{T'^)) = L^{M*) ® L^iT'^) . 
U is an isometry and U^^ given through 

(^Y"V)(^)-/ dfce'"V(fc,N) (11) 

JM* 

satisfies U^^Uip — ip for ip e iS(M'*), as can be checked by direct calculation. 
U^^ extends to an isometry from Tir to i^(R''). Hence must be injective 
and as a consequence U must be surjective, thus unitary. 

In order to determine the Zak transform of operators like the full Hamil- 
tonian in we need to discuss how differential and multiplication operators 
behave under the Zak transform, see |Bli| , |Zaj . Let P = —iVx with domain 
H'^(R.'^) and Q muhiplication by X on the maximal domain. Then 

UPU-^ = l®-iVP°'^ + fc0l, (12) 
UQU-^ = iVI, (13) 

where — iVP"^' is equipped with periodic boundary conditions or, equivalently, 
operating on the domain H'^{T'^). The domain of iV^ is Hr r\ H^^^iW^ , L'^iT'^)), 
i.e. it consists of distributions in H^{M* , L'^iT'^ j) which satisfy the y-dependent 
boundary condition associated with @ . In addition to H12|) and we notice 
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that multiplication with a F-periodic function like Vr is mapped into multipli- 
cation with the same function, i.e. L{Vr{x)U~^ = 1 ^Vr{y)- 

For later use we remark that the following relations can be checked using 
the definitions Q and ifTT)) . 

Remark 1. The Bloch-Floquet transform is usually defined as 

(t/i/')(fc,a;) ^e-'^-V(a; + 7), {k,x)eR^'^. (14) 

for -0 G In contrast to |7J|, functions in the range ofU are periodic in k 

and quasi-periodic in y, 

(Uijj) {k, y + 7) = c"="' (WV') (k, y) for all 7 e F , (15) 

{Ui^){k + j*,y) = {l{tl;){k,y) for all 7* e F* . (16) 

Our choice of using the Zak transform U instead of U comes from the fact that 
the transform of the gradient has a domain which is independent of /c G M* , see 
H12|) . As we shall see, this is essential for the application of the pseudodifferential 
calculus of Appendix A. <C> 

For the Zak transform of the free Hamiltonian one finds 
UHp„U-^ = / dkHpc,{k) 

with 

^pcr(fc) - ^(-iV^ + fc)' + l^r(y), fceM*. (17) 

For fixed k G M* the operator Hp^,{k) acts on L'^{T'^) with domain ij2(T'*) 
independent of fc G A/*, whenever the following assumption on the potential is 
satisfied. 

Assumption Ai. We assume that Vy is infinitesimally bounded with respect 
to -A and that cj) G C^f (M'*, M) and Aj G C{f (R'^, R) for any j G {1, . . . , d}. 

Here C^{M.'^, R) denotes the space of bounded smooth functions with deriva- 
tives of any order bounded. From this assumption it follows in particular that 
also the full Hamiltonian i?^ of l(2Jl is self-adjoint on iJ^(IR'^). Assumption (Ai) 
excludes the case of globally constant electric and magnetic field. However, since 
we are not concerned with the spectral analysis of 7?'^, but with the dynamics 
of states for large but finite times, locally constant fields serve us as well. 

The band structure of the fibered spectrum of iJpcr is crucial for the follow- 
ing. The resolvent Rl = {Ho{k) - A)"i of the operator iJo(fc) = K - iVy + k)^ 
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is compact for fixed k e M*. Since, by assumption, i?AVr is bounded, also 
R\ = (i?pcr(fc) — A)"^ — R^+R\VrRl is compact. As a consequence ifpor(fc) has 
purely discrete spectrum with eigenvalues of finite multiplicity which accumulate 
at infinity. A more detailed discussion can be found e.g. in ,Wi- . For definiteness 
the eigenvalues are enumerated increasingly as Ei{k) < E2{k) < i?3(fc) < . . . 
and repeated according to their multiplicity. The corresponding normalized 
eigenfunctions {ipn{k)}n&fi C H'^{T'^) are called Bloch functions and form, for 
any fixed k, an orthonormal basis of L^(T''). We will call En{k) the n}^ band 
function or just the n}^ band. Notice that, with this choice of the labelling, 
Enik) and ipn{k) are generally not smooth functions of k due to eigenvalue 
crossings. Since 

Hp.r{k - 7*) = r(7*) i/pe,(fc) r(7*)-i , (18) 
the band functions En{k) are periodic with respect to T* . 

Definition 2. A family of Bloch bands {i?„(fc)}„gi, X = [/_,/+] p| N, is called 
isolated, or satisfies the gap condition, if 

^M^dist( U{£;„(fc)}, U [Eraik)}) Cg > 0. 

In the following we fix an index set I C N for an isolated family of bands. 
Let Pi{k) be the spectral projector of Hpci {k) corresponding to the eigenvalues 
{En{k)}nGX, then Pj := ^^.j, dk Pi{k) defines the projector on the given isolated 
family of bands. In terms of Bloch functions Pi{k) = '^ni^i\'^n{k)){'^n{k)\. 
However, in general, ipn{k) are not smooth functions of k at eigenvalue crossings, 
while Px{k) is a smooth function of k because of the gap condition. Moreover, 
from H18|) it follows that 

Pi{k-j*) = rh*)Pi{k)T{r)-'- 
For the mapping to the reference space we will need the following assumption. 

Assumption A2. // the isolated family of bands {-En(^)}nei ^■s degenerate, 
in the sense that £ = |X| > 1, then we assume that there exists an orthonor- 
mal basis o/RanPi(fc) whose elements are smooth and T-equivariant 
with respect to k, i.e. '0j(fc — 7*) = ''"(7*)'0j(fc) for all j G {1, andj* G F*. 

In the case of a single isolated i-iold degenerate Bloch band (i.e. i?„(fc) = 
E^{k) for every rt G Z, \I\ — £), Assumption (A2) is equivalent to the exis- 
tence of an orthonormal basis consisting of smooth and r-equivariant Bloch 
functions. On the other side, if there are eigenvalue crossings inside the family 
of bands. Assumption (A2) requires only that tpji^) is an eigenfunction of the 
corresponding eigenprojection Pi{k) and not of the free Hamiltonian ffper(fc). 

From the geometrical viewpoint Assumption (A2) is equivalent to the trivi- 
ality of a complex vector bundle over T'', namely the bundle of the null spaces 
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of 1 — Pi{k) for k £ M*. In this geometrical perspective it is not difficult to 
see that Assumption (A2) is always satisfied if either d — 1 or £ — I. Indeed, 
classification theory for bundles implies that any complex vector bundle over 

= is trivial. As for ^ = 1, it is a classical result, due to Kostant and 
Weil, that smooth complex line bundles are completely classified by their first 
integer Chern class. In our case, the time-reversal symmetry of iJper implies 
the vanishing of the first integer Chern class, and therefore the triviality of the 
bundle. The same, and indeed slightly stronger, results can be proved with 
analytical techniques, as in |Ne| and references therein. By pushing forward 
the geometrical approach above, we expect that Assumption (A2) is generically 
satisfied for d < 3, as it will be discussed in |Pa| . 

In the presence of a strong external magnetic field the Bloch bands split 
into magnetic sub-bands. Generically, their first Chern number does not vanish 
and therefore Assumption (A2) fails. As well understood and discussed in the 
introduction, the nonvanishing of the first Chern number is directly linked to 
the integer quantum Hall effect |TKNNI ISi] . hence our interest in extending 
Theorem 3 to magnetic Bloch bands. The required modifications of our theory 
will be discussed in IPST3I . 



3 Space-adiabatic perturbation for Bloch bands 

Let Pn{k) = |iy9„(fc))(iy9„(fc)|. Then the projector on the n}^ band subspace is 
given through P„ = dk Pn{k). By construction the band subspaces are 
invariant under the dynamics generated by iJpor, 



^ 1 -'^ n ^ 



iE„{k)s^ P„ = for aU n e N , s e' 



Notice that P„ is not a spectral projector of i?pcr, in general, since in more than 
one space dimension it can happen that e.g. -E„(fc) < En+i{k) for all k e M* 
but inffc£'„+i(fc) < supj.£'„(fc). According to the identity H12(l . in the original 
representation i/por acts on the n"' band subspace as 

iJpcr^ = U-\En{k) (E>1)U^ = En{-iV^) , 

where %l} £ U~^Pn U L'^{W^). In other words, under the time evolution generated 
by the periodic Hamiltonian wave functions in the n**^ band subspace propa- 
gate freely but with a modified dispersion relation given through the n'^ band 
function En{k). 

In the presence of non-periodic external fields the subspaces P„7ir are no 
longer invariant, since the external fields induce transitions between different 
band subspaces. If the potentials are varying slowly, these transitions are small 
and one expects that there still exist almost invariant subspaces associated with 
isolated Bloch bands. To construct them, and to study the dynamics inside 
these almost invariant subspaces, we apply adiabatic perturbation to perturbed 
Bloch bands. 



9 



We first present a theorem which summarizes the main results of this sec- 
tion. The remaining parts give the results and the proofs of the three main 
steps in space-adiabatic perturbation theory: In Section 3.1 we construct the 
almost invariant subspaces associated with isolated Bloch bands. In Section 3.2 
we explain how to unitarily map the decoupled subspace to a suitable refer- 
ence Hilbert space. In this reference space the action of the full Hamiltonian 
is given through a semiclassical pseudodifferential operator, whose power series 
expansion can be computed to any order in e. This effective Hamiltonian is con- 
structed in Section 3.3 and we compute explicitly its principal and subprincipal 
symbol. The main technical innovation necessary in order to apply the scheme 
to the present case is the development of a pseudodifferential calculus for op- 
erators acting on sections of a bundle over the flat torus M* , or, equivalently, 
acting on the space Ht- This task is deferred to Appendix A. 

Before going into the details of the construction we present a theorem which 
encompasses the main results of this section. Generalizing from (|10|l it is con- 
venient to introduce the following notation. For any separable Hilbert space Tif 
and any unitary representation r : F* — > l4{Ti.f), one defines the Hilbert space 

equipped with the inner product 

^>L?= / dfc(V;(fc), ^(fc))w,. 
JM' 

Using the results of the previous section and imposing Assumption (Ai), the 
Zak transform of the full Hamiltonian in ([SJ is given through 

Hi ■.= UH'U-^ = ]^i^^iVy + k-A(ieVl))\vr{y)+^{}eVl) (19) 

with domain Ll{W^ , {T'^)) . 

The application of space-adiabatic perturbation theory to an isolated family 
of bands {ii^n(^)}nei yields the following result, where the reference Hilbert 
space for the effective dynamics is JC := L'^{M*) ® with £ := dimPi(/c). 

Theorem 3 (Peierls substitution and higher order corrections). Let 

{En}nei ^6 isolated family of bands, see Definition\^ and let the Assump- 
tions (Ai) and (A2) be satisfied. Then there exist 

(i) an orthogonal projection H^ G B{Tir), 

(ii) a unitary map G B{Il^T-lr,IC), and 

(iii) a self-adjoint operator h G B{IC) 
such that 

II |hO(e°°), II - Pi II = 0(e) 
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and 

rH^t _ jje* ^-M jje^^j^e || ^ 0{e°° {I + \t\)) . 



e 



The effective Hamiltonian h is the Weyl quantization of a semiclassical symbol 
h € S^^^{e, B{C^)) with an asymptotic expansion to any order. The B{C^)- 
valued principal symbol h^{k^r) has matrix- elements 

hoik,r)a,p = {Mk- Air)), Hoik, r)Mk-A{r))), (20) 

where a, P € {1, ■■■,(} and H^ik, r) is defined in \24^ . 

The general formula for the subprincipal symbol of the effective Hamiltonian 
can be found in |PSTi| . The structure and the interpretation of the effective 
Hamiltonian are most transparent for the case of a single isolated band. 

Corollary 4. For an isolated i-fold degenerate eigenvalue the B(C^)- 

valued symbol ft,(fc,r) ~ hoik,r) + ehiik,r) + Ooie'^) constructed in Theorem\^ 
has matrix- elements 

hoik, r)^p = {Eik - Air)) + 0(r))(5,^ (21) 

and 

hiik,r)^fi = -{- V(^(r) + V^(fc) x B(r)) • Aik)^p - Bir) ■ M(fc),0 (22) 



a,0(r) - diEik) (d.Aiir) - diA,ir))j ^■(fc)„^ 
~{d,Ai - diA,)ir) Re [i (di^^ik), (i/p„. - E)ik) d,^p(k))^^ 

where summation over indices appearing twice is implicit, kik,r) = k — Air), 
and a, (3 G {1, . . . ,f\ . The coefficients of the Berry connection are 

Aik)^p^\{,p^ik),Vi;pik))^^. (23) 

In dimension d — i the subprincipal symbol H22|l has a straightforward phys- 
ical interpretation. The 2-forms B and M are naturally identified with the 
vectors B = curM and 

Mik)^p = i (VVa(fc), x(iJp„.(fc) - Eik))Vijpik))^^ . 

Therefore the symbol of the effective Hamiltonian has the same form as the 
energy of a classical charge distribution in weak external fields, in first order 
multipole expansion. In this sense ■4(fc) is interpreted as an effective electric 
dipolc moment and A^(fc) as an effective magnetic dipole moment. 

Remark 5. Our results hold for arbitrary dimension d. However, to simplify 
presentation, we use a notation motivated by the vector product and the duality 
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between l-forms and 2-forms for d = 3. If d 7^ 3, then B, fin and M„ are 2- 
forms. The inner product of 2-forms is 

d d 

where * denotes the Hodge duahty induced by the eucHdian metric, and for a 
vector field w and a 2-form F the "vector product" is 

d 

iw X F), {*-\w A *F))j = ^ u-,F,, , 

i=l 

where the duahty between l-forms and vector fields was used implicitly. <C> 

Theorem O is a direct consequence of the results proved in Propositions El El 
andjHl The proof of CoroUary^lis given at the end of this section. 

As mentioned before, the main idea of the proof is to adapt the general 
scheme of space-adiabatic perturbation theory to the case of the Bloch electron. 
While formally this seems straightforward, one must overcome two mathemat- 
ical problems. First of all, in the present case the symbols are unbounded- 
operator- valued functions. One can deal with unbounded-operator- valued sym- 
bols by considering them as bounded operators from their domain equipped 
with the graph norm into the Hilbert space, see e.g. |DiSj| . The second, more 
serious problem consists in setting up a Weyl calculus for operators acting on 
spaces like i^(R'^,7if). This is done in Appendix A and we will use in this 
section the terminology and notations introduced there. 

The results of Appendix A allow us to write the Hamiltonian iJ| as the Weyl 
quantization HQ{k,iey k) of the r-equivariant symbol 

Ho{k,r) = ^{-iV^ + k-A{r)y + Vr{x)+cb{r) (24) 

acting on the Hilbert space Hi L^(T'*,dx) with constant domain V := 
iJ^(T'^). For sake of clarity, we spend two more words on this point. For 
any fixed {k, r) € K^'', -ffo(fcj t) is regarded as a bounded operator from V to Hi 
which is r-equivariant with respect to the bounded representation ti := t|xj act- 
ing on D and the unitary representation T2 := r acting on Tif , see Definition 12 II 
Then the general theory developed in Appendix A can be applied. The usual 
Weyl quantization of Hq is an operator from S'{M.'^,T>) to S'{W'-,Hi) given by 

i?o = i( - iVj, + fc - A(i£Vfc))' + Vviv) + (^(ieVfe) . (25) 

Then Hq can be restricted to Lf^^[M.'^,V), since A and (p are smooth and 
bounded. Since is a r-equivariant symbol, Hq preserves r-equivariance and 
can then be restricted to an operator from L^(M'^,X') to L'^{W^,Hi). To con- 
clude that ij^ . restricted to L^iW'-.V), agrees with ((T^ . it is enough to recall 
that iV^ is defined as iVfc restricted to H^OHr and to use the spectral calculus. 
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Moreover, if one introduces the order function w{k,r) (1 + fc^), then 
Hq G S'!^{B{'D,Ti.)). More generally, we will give the proofs for any symbol 
H G -S'^(e, B{V, H)), whose principal symbol is then denoted by -ffo- 

3.1 The almost invariant subspace 

In this section we construct the adiabatically decoupled subspace associated 
with an isolated Bloch band. Similar constructions have a considerable history 
and we refer to [MaSoL [ReSoL |PSTi| |Tei| and references therein. 

Given an isolated family of bands {En{k)}n£x, we define 7ro(fc, r) = Px{k — 
A(r)). It follows from the r-equivariance of Hq and from the gap condition 
that TTo e S^{B{ni)). We also define the shorthand A{e) = Oo(e"), where the 
subscript expresses that a family A{e) G B{H) is 0(e") in the norm of bounded 
operators. By A{e) = Oo{£°°) we mean that A{e) = C'o(e") for any n G N. The 
remaining notation is defined in Appendix A. 

Proposition 6. Let {En}nei be an isolated family of bands and let Assumption 
(Ai) be satisfied. Then there exists an orthogonal projection T¥ G B{T-Lr) such 
that 

=Oo(£°°) (26) 

and II"^ = TT + 0{e°^), where tt is the Weyl quantization of a r-equivariant 
semiclassical symbol 

^x^e^vr, in Sl{e, B{ni)) , 

whose principal part 7ro(fc,r) is the spectral projector of HQ{k,r) corresponding 
to the given isolated family of bands. 

Proof. Wc first construct tt on a formal symbol level. 

Lemma 7. Let w{k,r) = (1 + fc^). There exists a unique formal symbol 

oo 

TT = ^ e^'ttj- g Ml{e, B{Hi)) D M^{e, B{Hi, V)) 

3=0 

such that 7ro(fc, r) — Px{k — A{r)^ and 

(i) TrjjTT = TT, 
(ii) TT* = TT, 

(iii) i?tt7r-7rtJi7 = 0. 

Proof. We construct the formal symbol vr locally in phase space and obtain by 
uniqueness, which can be proved as in PSTi , a globally defined formal symbol. 

Fix a point zq = (fco,ro) G R^''. From the continuity of the map z ^ H{z) 
and the gap condition it follows that there exists a neighborhood Uza of zq 
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such that for every z G Wz^, the set {En{z)\n^x can be enclosed by a positively- 
oriented circle t^izo) C C independent of z in such a way that A(zo) is symmetric 
with respect to the real axis, 

dist(A(zo),cr(iJ(z))) > for all z C^U^^ (27) 

and 

Radius(A(zo)) < ■ (28) 

The constant Cg appearing in (|27|l is the same as in Definition [21 and the exis- 
tence of a constant Ci independent of zq such that (|28(l is satisfied follows from 
the periodicity of {En{z)}nei and the fact that A and (f) are bounded. Indeed, 
A can be chosen r*-periodic, i.e. such that A(fco -|- 7*,7'o) — A(A:o,ro) for all 

7* € r*. 

Let us choose any C e A(zo) and restrict all the following expressions to 
z € UzQ- We will construct a formal symbol R{C) with values in B{Tii^'D) — 
the local Moyal resolvent of — such that 

{H~0^R[0^1n, and i?(C) (i? - C) = li? on . (29) 
To this end let 

where according to (|?7jl Rq{Q{z) G B[Hi,'D) for all z E Uzq, and, using differ- 
entiability of H{z), d"Ro{C){z) e B{Hi,'D) for all z e U^g. By construction one 
has 

(i/-C)«i?o(C)-lWf +Oo(e), 

where the remainder is 0{e) in the ;B(7if)-norm. We proceed by induction. 
Suppose that 

n 

R(-){C) = Y,e^R,{0 
j=o 

with Rj{(){z) e B{H{, T>) for all z E Uzg satisfies the first equality in (|29|l up to 
©(£"+1), i.e. 

(i/ - C) tt i?("HC) = Iw, + e"+'i?„+i(C) + a(£"+') , (30) 
where En+i{C){z) E B{Hf). By choosing 

i?„+l(C)--i?o(C)£^n+l (31) 

we obtain that i?("+i)(C) = i^^'^HC) + e"+^i?„+i(C) takes values in B{nf,V) 
and satisfies the first equality in 129|l up to ©(e"''"^). Hence the formal symbol 
R(C,) = J^JLo -^j (0 constructed that way satisfies the first equality in H29|l 
exactly. By the same argument one shows that there exists a formal symbol 
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R{() with values in B{Hi, V) which exactly satisfies the second equality in (j^ . 
By the associativity of the Moyal product, they must agree: 

i?(C) = i?(C)tt(^-C)tt^(C) = i?(C) onu.^. 

Equations H29|) imply that i?(C) satisfies the resolvent equation 

i?(C)"i?(C') = (C-C')^(C)tti?(C') onW,„ (32) 

for any C,X' G A(zo)- From the resolvent equation it follows as in that 
the S(Ht, 2?)- valued formal symbol vr ~ '^JLo^^'^j defined through 

{z)-=^ i dC Rj (C, z) on (33) 

satisfies (i) and (ii) of Lemma [7| As for (iii) a little bit of care is required. Let 
J : V ^ Tif he the continuous injection of T> into Tif . Using H33|) and (|32|l it 
follows that TT JjJi?(C) = i?(C)^tt7r for all C, € A(zo). Moyal- multiplying from 
left and from the right with H — C, one finds H J = Jtt'^H as operators in 
B{'D, H[). However, by construction H ^tt takes values in B{7i{) and, by density 
of V, the same must be true for tt'^H. 

We are left to show that tt e M^ie, S(Hf)) n M^'{e, S(Hf , V)). To this end 
notice that by construction tt inherits the r-equivariance of H, i.e. 

nj{k - j*,q) = r(7*)7rj(fc,g)T(7*)-i . 

From H33|) and (|28|l we conclude that for each a G N^'' and j G N one has 

||(9>,)(z)|| <2^a. sup ||(a,"i?,)(C,z)l|, (34) 
CeA(zo) 

where || • |j stands either for the norm of B{H{) or for the norm of B{Ti,i, V). In 
order to show that tt G M^(e, B{H[)) it suffices to consider z — (fc, r) G M* x 
since t(7*) is unitary and thus the ;B(7if)-norm of tt is periodic. According to 
1)34(1 we must show that 

||(a,"i?,)(C,z)|i5(>,,) < Vz G C e A(zo) (35) 

with Caj independent of zq G M* x R''. 

We prove (|35|l by induction. Assume, by induction hypothesis, that for any 
j < n one has that 

i?,(C) e 5^(B(7^f)) n5r(S(Wf,25)) (36) 

uniformly in (, in the sense that the Frechet semi-norms are bounded by C- 
independent constants. Then, according to Proposition]^ En+i{(), as defined 
by l|5n|) . belongs to (S(7if)) uniformly in By r-equivariance, the norm of 
En+i{C) is periodic and one concludes that £^„+i(C) G Sl{B{Hf)) uniformly in 
(. It follows from (|31|l that H36|) is satisfied for j = n + 1. 
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We are left to show that H3t)|l is fulfiUed for j — 0. We notice that according 
to (EH) one has for aU z e M^'* 

\\RoiO\Wn, - miz) - CriWn, - dist(C, ^ ^ " 

By the chain rule, 

||(a,i?o)(C,^)lb(w,) = ||(i?o(C)(a.i/o)i?o(C))WllB(«,)- (37) 

Since OzHq -Ro(C) is a r-equivariant i3(7if)-valued symbol, its norm is periodic. 
Therefore it suffices to estimate its norm for z £ M* x K'*, which yields the 
required bound. For a general a G N^'', the norm of 9"i?o(C) can be bounded 
in a similar way. This proves that i?o(C) belongs to S].{B{7if)) uniformly in C^. 
On the other hand 

\\Ro{k,r)\\Bin„v) = I|(l + A,)i?o([fc]+7*,r)||H(w,) 

= ||(l + A,)r(7*)i?o([fc],r)r-i(7*)||5(«,) 

< C7||(l+7*')(l + Ax)i?o([fc],r)||B(«,) 

< C"(l+7*^) < 2C"(l + fc2), 

where we used the fact that ||(1 + A2,)i?o(^) HsCHf) is bounded for z e M* x M.'^. 
The previous estimate and the fact that OzHq i?o(C) £ S].{B{Tii)) yield 

l|(a,i?o)(C,^)ll \\{RoiC)idzHo)Ro{0) 
< C(l + fc2). 

Higher order derivatives, are bounded by the same argument, yielding that 
i?o(C) belongs to S^{B{'Hf,'D)) uniformly in (. This concludes the induction 
argument. 

From the previous argument it follows moreover that 

||(5,"i?,)(C, z)\\b(h,.v) < w{z) yzeUz„,C(^ A(zo) (38) 



and concludes the proof. □ 



with Ccj independent of zq G R^''. By this implies tt G MJf {£,B{7iu'D)) 



Proof of Proposition From the projector constructed in Lemma H one 
obtains, by resummation, a semiclassical symbol tt G S].{e,'Hi) whose asymp- 
totic expansion is given by X]j>o ^"'"'i- Then according to Proposition 1251 Wevl 
quantization yields a bounded operator tt G B{l-Lr), which is approximately a 
projector in the sense that 

=^? + Oo(e°°) and ^* 

We notice that Proposition El implies that hIt: G 5^' (e, BCWt))- But t- 
equivariance implies that the norm is periodic and then ft tt belongs indeed 
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to Sl{e,B{n{)). Then tt^H = {H^tt)* belongs to the same class, so that 
[ff, 7r]j G Sl{£,B{'Hf)). This a priori information on the symbol class, together 
with Lemma [7|(iii), assures that 

[i/,^] =Oo(£°°) (39) 

with the remainder bounded in the yS(7iT-)-norm. 

In order to get a true projector, we proceed as in |NeSo| . For e small enough, 

let 

n^:=f/ dC(5f-C)-^ (40) 

Then it follows that W"^ ,W = ^ + Oa{£°°) and 

II [i?, ] < C\\ [H, ^ ||e(„o = O(e-) . 

□ 



3.2 The intertwining unitaries 

After having determined the decoupled subspace associated with an isolated 
family of Bloch bands, we aim at an effective description of the intraband dy- 
namics, i.e. the dynamics inside this subspace. In order to get a workable 
formulation of the effective dynamics, it is convenient to map the decoupled 
subspace to a simpler reference space. The natural reference Hilbert space for 
the effective dynamics is /C := L'^[T'^*)(Si<C\ where £ := dimPx(fc) and T''* is M* 
with periodic boundary conditions. Notation will be simpler in the following, 
if we think of the fibre as a subspace of Tif. In order to construct such a 
unitary mapping, we reformulate Assumption (A2). 

Assumption A.^. Let {-E'n(A:)}„gx be an isolated family of hands and let tTj. G 
B{Hf) be an orthogonal projector with dimTTj. = £. There is a unitary- operator- 
valued map uq : M^'' —> 14 (Tif) so that 

uo{k,r)TTo{k,r)ul{k,r) ^ TT^ (41) 

for any {k, r) G M^'', 

wo(fc + 7*,r) = Mo(fc,r)r(7*)-i , (42) 
and Uq belongs to S^{B{Ti.{)). 

Clearly, 

^o(fc + 7*,'') = t(7*)wo(^>0- (43) 

An operator- valued symbol satisfying H43(l (resp. 142|l ') is called left r-covariant 
(resp. right r-covariant). 

The equivalence of (A2) and (A^ can be seen as follows. According to 
Assumption (A2), there exists an orthonormal basis {V'j(fc)}j=i of RanPi(fc) 
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which is smooth and r-equivariant with respect to k. Let tTi- 7ro(fco,ro) for 
any fixed point (fco,ro). By the gap condition, diniTTr — diniPi(fc). Then for 
any orthonormal basis {Xj}j=i for RanTTr, the formula 

e 

Mk,r) :=5]|x,)(^,(fc-A(r))| (44) 
i=i 

defines a partial isometry which can be extended to a unitary operator wo(fc, r) G 
U{Ti.{). The fact that spans RanPi(fc) implies (|41|l . and the r- 

equivariance oi'4'j{k) reflects in (02) • 

Viceversa, given mq fulfilling Assumption (A2), one can check that the for- 
mula 

■il)j{k ~ A{r)) ■.= u*Q{k,r)xj, 

with {xjYj^i spanning RanTTj., defines an orthonormal basis for RanPi(fc) which 
satisfies Assumption (A2). 

After these remarks recall that the goal of this section is to construct a 
unitary operator which allow us to map the intraband dynamics from Rann"^ 
to an £-independent reference space JC C Tircf- Since all the twisting of Tir 
has been absorbed in the r-equivariant basis {V'jlj^ii or equivalently in uq, the 
space Tircf can be chosen to be a space of periodic vector-valued functions, i.e. 

We introduce the orthogonal projector Ifr := vTr G B{Hicf) since the effective 
intraband dynamics can be described in 

/C := Rann,. ^ L^=i(M'', = L^(T^*, 

as it will become apparent later on. Recall that £ = dimPi(fc) = dimTTf. 

Proposition 8. Let {En}nei be an isolated family of bands and let Assumptions 
(Ai) and (A2) be satisfied. Then there exists a unitary operator : Tir — * '^ref 
such that 

jje n-^ * = n^. (45) 

and — it + OQ{e°°), where u x X]j>o^"'"i belong to S^{e,B{H{)), is right 
T-covariant at any order and has principal symbol uq. 

Proof. By using the same method as in Lemma 3.3 in |PSTi| , one constructs 
first the formal symbol X]j>o^"'"i- Since uq is right r-covariant, one proves 
by induction that the same holds true for any Uj. Indeed, by referring to the 
notation in |PSTi| , one has that 

Un+l = (fln+l + bn+l)uo 

with a„+i = —^An^i and fen+i = [tt,., Bn+i]. From the defining equation 
y(") _ 1 ^ + C'(e"+2) 
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and the induction hypothesis, it fohows that An+i is a periodic symbol. Then 
:— + £"+ia„+iMo is right r-covariant. Then the defining equation 

I TT I - TT, = e"+iB„+i + C'(e"+2) 

shows that Bn+i is a periodic symbol, and so is fo„+i. Hence is right t- 
covariant, and there exists a semiclassical symbol u x '^j£-'uj so that m g 
5i(e,S(Hf)). 

One notices that right r-covariance is nothing but a special case of (ti,T2)- 
equivariance, for T2 = 1 and ri = r. Thus it follows from Proposition 1251 that 
the Weyl quantization of u is a bounded operator u G B{Ht, Tiief) such that: 

(i) uu* ^ lH„f + Co(e°°) and u*u = 1^^ + Oo(e°°), 

(ii) Mffu* =ni. + C'o(e°°). 

Finally we modify u as in IPSTil by an Oo(£°°)-term in order to get the 
unitary operator If^ G U{Hr,'Href)- □ 

3.3 The effective Hamiltonian 

The final step in spacc-adiabatic perturbation theory is to define and compute 
the effective Hamiltonian for the intraband dynamics and to compute its lower 
order terms. This is done, in principle, by projecting the full Hamiltonian Hy^ 
to the decoupled subspace and afterwards rotating to the reference space. 

Proposition 9. Let {-Enjnei be an isolated family of bands and let Assumptions 
(Ai) and {A2) be satisfied. Let h be a resummation in S^^j^{e, B {?{[)) of the 
formal symbol 

h = uiTriHi7T^u* e Ml^^{e,B(Hi)). (46) 
Then h G S(7^rof), \h,^A ^ and 

(e^i^z* _ e"'''* U')W ^ 0o(e°°(l + \t\)) . (47) 

Remark 10. The definition of the effective Hamiltonian is not entirely unique 
in the sense that any H^ff satisfying H47I) would serve as well as an effective 
Hamiltonian. However, the asymptotic expansion of Hcs is unique and therefore 
it is most convenient to define the effective Hamiltonian through H4()|l . <> 

Proof. In the proof we denote iJgp as H to emphasize the fact that it is the 
Weyl quantization of i? G ^^(e, B{V, Hi)). 

First note that (|46|) follows from the following facts: according to Lemma[7| 
and Proposition 1241 we have that 

^ ft TT G Mf{e, BiHf)) = m;(£, BiUf)) , 

where we used that r is a unitary representation. With Proposition|Slit follows 
that h G M^^^{e,B{ni)). Therefore h G -B(Wrcf) follows from Proposition!^ 
while [h, Hi ] = is satisfied by construction. 
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It remains to check (I47II : 

(e-'^*-i7"e-'''*C/^)ff = (e^'^*-e-'^'*''^'*)?f + Oo(e°°) 

^ (e~-^-*-e-'^'*^^'*)7? + Oo(£°°) 

where the last equahty foUows from the usual Duhammel argument and the fact 
that the difference of the generators is Oo{e°°) in the norm of bounded operators 
by construction. □ 

Since [/ijllr] = 0, the effective Hamiltonian will be regarded, without dis- 
tinctions in notation, either as an element of S(7iref) or as an element of B{IC). 

We compute the principal and the subprincipal symbol of h for the special 
but most relevant case of an isolated eigenvalue, eventually ^-fold degenerate, 
i.e. En{k) = E{k) for every n £ I, \2\ — i. Recall that in this special case 
Assumption {A2) is equivalent to the existence of an orthonormal system of 
smooth and r-equivariant Bloch functions corresponding to the eigenvalue E{k). 
li i = 1 then Assumption (A2) is always satisfied. The part of uq intertwining 
TTg and TTi- is given by equation H44|l where ^j{k) are now Bloch functions, i.e. 
eigenvectors of Hpc-[{k) with eigenvalue E{k). 

Proof of Corollary^ In the following h is identified with TTj^hnj- and regarded 
as a ,B(C'')-valued symbol. We consider the matrix elements 

h{k,r)af3 ■■= {Xa,h{k,r)x0) 

for a, f3 G {1, . . . , ^}, where we recall that Xa ~ UQ^k, r)'ijja {k — Mr)). Equation 
H21|) follows immediately from the fact that /iq = uq ^0 w|5 and that i/iq are 
Bloch functions. As for /ii, we use the general formula of |PSTi| , which reads, 
transcribed to the present setting, as 

/jia/3(A:,r) = -i(i^„(:fc), {^(fc) + 0(r),?/>/3(fc)}) 

- i(^„(fc),{(i/p„.(fc) - E(k)),i,fi(k)}) . (48) 

Here {A, Lp] ~ A • Vfe(p — VfcA • Vr'p are the Poisson brackets for an operator- 
valued function A(k^r) acting on a vector- valued function ip(k,r). We need to 
evaluate l|48|l . Inserting H44|l and performing a straightforward computation the 
first term in (|48|) gives the first term in (|22|l while the second term contributes 
to the a(3 matrix element with 

d 

J2 {djAi - diA,){r) {Mk),di{Hp,r - E){k) d,Mk))n, ■ 
3,1=1 

The derivative on (iJpcr — E) can be moved to the first argument of the inner 
product by noticing that 

= V(Va, (-ffpcr - E)<j,) = ( VV'a, (Fpor " E)q)) + V(i7per " E)q^) 
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since ipa is in the kernel of (i^pcr — E). Finally the imaginary part of 

d 

^ J2 (djAi - diA,){r) (diMk), (i^per - E)Ck) d.M^H, 

3,1=1 

vanishes, as can be seen by direct computation, concluding the proof. □ 



4 Semiclassical dynamics for Bloch electrons 

We have now at our disposal the tools to establish the link between the Schrodin- 
ger equation |2Jl and the corrected semiclassical equations of motion © • To this 
end we specialize to the case of a non-degenerate Bloch band £"„. The phase 
space for (|SJ) is R'' x R'', since we use the extended zone scheme, and we denote 
by $* the corresponding solution flow. Since the effective Hamiltonian is written 
in canonical variables, it is necessary to switch in l^l to (r, k) with k — k + A{r). 
In the new coordinates the solution flow is denoted by '^^ and 

¥*(r,A;) = ($*,(r,A;- A(r)), ,(r, fc - A(r)) + A(r)) . 

Let us consider any admissible semiclassical observable a = a(ex, — iV^;) acting 
on the "physical" Hilbert space L^(R^,dx). Its symbol is transported by '^^ 

to a o $* with Weyl quantization a o . On the other hand the operator a is 
transported by the Heiscnberg equation as e'^^*/^ae~'''^^*/^. Our assertion is 
that on the subspace W^L'^{W^) := U^*WHr, 11= and as constructed in the 
previous section, these two operators are uniformly close to order . 

Theorem 11. Let En be an isolated, non- degenerate Bloch band, see Defini- 
tion\^ and let the potentials satisfy Assumption (Ai). Let a S C^(R^'') be 
T* -periodic in the second argument, i.e. a(r, fc + 7*) — a(r,k) for all 7* S T* , 
and a = a(sx, — iVa;) be its Weyl quantization. Then for each finite time-interval 
/ C M there is a constant C < 00 .such that for t £ L 

In particular, for ^jq S II^W we have that 

I (Vo, e'«^*/= a e-'^'*/= V'o) - (^0, V'o) \ < C Uof ■ (49) 

Theorem llll is an Egorov-type theorem, see |Ro| . An unconventional feature 
is that the flrst order corrections are treated by considering an e-dependent 
Hamiltonian flow instead of having a separate dynamics for the subprincipal 
symbol of an observable. 

By exploiting the relation between Weyl-quantized operators and Wigner 
transforms, one can easily translate H49(l to the language of Wigner functions. 
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For a detailed discussion on how Theorem relates to alternative approaches 
to the semiclassical limit in perturbed periodic potentials we refer the reader to 

ma- 

To prove Theorem our strategy is to first establish a corresponding 
Egorov theorem in the reference space and then to pull back to L^(K'^,dx). 

Proposition 12. Let E he an isolated non- degenerate Block hand and let h be 
the effective Hamiltonian constructed in Theorem\^ which acts on the reference 
space JC Ll^^iR'^) of T* -periodic L'^^^- functions. Let $* : M^rf ^ ^2d ^/^g 
Hamiltonian flow generated hy the Hamiltonian function 

hc\(k,r) = hQ{k,r) +ehi{k,r) . 

Then for any semiclassical ohservable a ~ ao(A:,ieVfe) + £ai(fc,i£Vfe) with a G 
S^{e,C) we have that 



a o$t II <CTe^ 
uniformly for any finite interval in time [—T, T] . 



(50) 



Proof. Since the Hamiltonian function is bounded with bounded derivatives, it 
follows immediately that a o cf)* g S^{e) and that ^ (a o $*) g S^{e). Therefore 
the proof is just the standard computation 



0$* 



dt' 



di' 

/i, (a o 



- -iL(ao$*-') 



-\ht je 



together with the fact that the integrand is Oie^) in the norm of bounded 
operators, since by construction 

and, computing the expansion of the Moyal product. 



/i, a o $* 



□ 



In order to obtain the Egorov theorem for the physical observables, we need 
to undo the transform to the reference space and the Zak transform. We start 
with the simpler observation on how the Zak transform maps semiclassical ob- 
servables. 
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Proposition 13. Let a 6 S^{e, C) be T* -periodic, i.e. a{r, k + 7*) — a{r, k) for 
all 7* G r*. Let b{k, r) ~ a{r, k) then b 6 5'^(£, C) and 

a = U*bU, 

where the Weyl quantization is in the sense of a ^ a{ex, — iVx) acting on i^(R'^) 
and b = b{k, eiS/k) acting on TCt. 

Remark 14. An analogous statement cannot be true for general operator- 
valued r-equivariant symbols. For example, the symbol b{k,r) := Hpci{k — 
A(r)) is r-equivariant and in particular a semiclassical observable. However, 
the corresponding operator in the original representation is 

WbU = -i( - iV:, - A{ex)f + Vr{x) 
which cannot be written as a e-pseudodifferential operator with scalar symbol. 

Proof. We give the proof for a(-,fc) G 5(R'^). The general result follows from 
standard density arguments, see |DiSj| . For tp G iS(M'*) we have according to 
H64(l the explicit formula 

(a(ex, -iV,)^) (x) - —1^ J2 I H^"") ^) e'^f^-^^/'e-"- V(x + 7) • 

(51) 

On the other hand for {l/(ip){k,r) ~: Lp{k,r) by definition it holds that 

(6(fc, ieVk)(p) (fc, r)^y2 [ dri{n) (7, ri) e-'^("-'^)/2e'''- " e^, r) . (52) 

The assumptions on a and ip guarantee that all the integrals and sums in the 
following expressions are absolutely convergent and thus that interchanges in 
the order of integration are justified by Fubini's theorem. 
We compute the inverse Zak transform of H52|l using Hll|l . 

{U-^bLp){x)^ (53) 
= V / dfc /" d77(J-6)(7,?7) e'*=-^e-'^(''-'^)/2e''^-V(fc-£'7, W) 

= V /" d?7(jf^6)(7,77)e'^(''-'^)/2e'"''-^ / dfce'('=-"'')-(^+''V(fc - £^7, N) . 

The r-equivariance of tp implies that the function f{k,y) :— e''^'^(/9(fc, [y]) is 
exactly periodic in the first variable. Then the integral in dfc can be shifted by 
an arbitrary amount, so that 

dfce'('=-^'')-(=^+^V('fc-e'7, N) = / dfce''=-(^+''V('fc, [a;-K7]) = VX2:-H7). 

M* Jm* 

Inserting this expression in the last line of (|53|1 and comparing with (|51|l con- 
cludes the proof. □ 
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Before we arrive at the proof of Theorem llll one has to study how the unitary 
map constructed in Section 3.2 maps observablcs in the Zak representation to 
observables in the reference representation. 

Proposition 15. Let b — 5o(fc,£iVfc) + £ fei(fc, eiVfc) with symbol b G S^{e,C) 
which is T* -periodic in the first argument. Let : H'^Ht fC be the unitary 
map constructed in Section 3.2. Then 

U'li^bW U"* =c + 0{e^), 

where c(e, fc, r) = {b o T) (fc, r) with 

T : M^d ^ R2d ^ ^(k + eA^m{k- A(r)) VA„.(r), r + eA{k - A(r))) . 

Here and in the following, summation over indices appearing twice is implicit. 

Proof. In order to compute c = u^ir'^b'^-K'^u* , observe that, since b is scalar- 
valued, the principal symbol remains unchanged, i.e. cq = wotto&ottoWq = 6o- 
For the subprincipal symbol we use the general transformation formula H48|l 
obtained for the Hamiltonian, which applies to all operators whose principal 
symbol commutes with ttq. In this case the eigenvalue E in (|48|l must be replaced 
by the corresponding principal symbol and a term for the subprincipal symbol 
bi must be added. Hence we find that 

ci(fc,r) = -i(^(fc-A(r)),{6o(fc,r),^(fc-A(r))}) 

+ {ijj{k - A(r)), bi{k, r)ij;{k - A{r))) 
= dkjoik, r) i (^(fc - A(r)), 9,„^(fc - A{r))) a„A„(r) 

+ drMk, r) i (^(fc - A(r)), 9„^(fc - A{r))) + b^(k, r) 
= dk,Mk,r)A„,{k- A{r))dnAn{r) 

+ drMk,r) An{k - A{r)) + bi{k,r) , 

where summation over indices appearing twice is implicit. Now a comparison 
with the Taylor expansion of (b o T^{k, r) in powers of e proves the claim. □ 

We have now all the ingredients needed for the 

Proof of Theorem ] 111 Let a S C^(K^'') be r*-periodic in the second argument, 
then according to Proposition^] we have 

n^e'^'*/^ae-'^'*/''n^ e'^^*/^6e-'^^*/^ff W (54) 

with b{k, r) ~ a{r, k). With Theorem 13 and Proposition^] we find that 

n^e'^^*/£5e-'^5*/"ff = C/"e'''*/"?e-"'*/^C/" + 0(£2), (55) 

where c{e,k,r) — (6oT)(fc,r). Now we can apply Proposition E] to conclude 
that 
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Since, for e sufficiently small, T is a diffeomorphism, one can write 

c o = c o T"^ o T o o T"^ oT =:co T^^ o$*or = 6o¥*oT, 

where the flow $* in the new coordinates will be computed explicitly below. 
Inserting the results into H55|) . one obtains 

jjsgiffjt/e^g-iffjt/sjje ^ (6o$* oT) + 0(£2) 

= (6o$*)ff + 0(e2), 

where we used Proposition E| for the second equality. Inserting into H54|l we 
finally find that 

e'^'*/" a e-'^'*/^ = (a o $*) + ^(e^) , (55) 

where we did not make the exchange of the order of the arguments in a explicit. 

Since the flow is determined only in approximation and only through its 
vector field, we make use of the following lemma. 

Lemma 16. Let : M^"* x M ^ K.^"* he the flow associated with the vector field 

Vi G C^(K2d^]R2d)^ i = 1,2. 

(i) If for all a G N^'' there is a Ca < 00 such that 

sup I d" {vi -V2)ix)\ < Ca , 

then for each bounded interval / C M there are constants C/.q < 00 such 
that 

sup \d''i'^{-^'2)ix)\<Ci,^e\ (57) 

tei,xeK^'' 

(ii) Let a G S^{e,C). If j5?| ) holds for the flows $1, <I>2, then there is a constant 
C < 00, such that for all t £ I 

||^°*!-^°^2|Ib(L^(M^)) <^£'- 

Proof. Assertion (i) is a simple application of Gronwall's lemma. Assertion (ii) 
follows from the fact that the norm of the quantization of a symbol in is 
bounded by a constant times the sup-norm of finitely many derivatives of the 
symbol, which are 0{e'^) according to (|57|l . □ 

According to assertion (ii) of the lemma it suffices to show that 

$'(r,fc)= ($*,(r,fc-A(r)), Jr, A: - A(r)) + A(r)) + ©(e^) 

in the above sense, where is the flow of Q. And from assertion (i) we infer 
that it suffices to prove the analogous properties on the level of the vector flelds. 
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Through a subsequent change of coordinates we aim at computing the vector 
field of up to an error of order 0{e^). We start with the vector field of 
The effective Hamiltonian on the reference space including first order terms 
reads 

h{r,k) = E{k - A{r)) + (l){r) (58) 
- £ (i^Lor(r, WE{k - A{r))) ■ A{k - A{r)) + B{r) ■ M{k - A(r))) , 
with the Lorentz force 

FLor(r, VE{k - A{r))) = -W<j){r) + VE{k - A{r)) x B{r) . 
Componentwise, the canonical equations of motion are 
r,- = dkMr,k) = dk,E{k - A{r)) 

-edk, (FLor(r, k - A{r)) ■ A{k - A{r)) + B{r) ■ M{k - A(r))) , 

kj = -dr,h{r,k) = -dj(l>{r) + diE{k - A{r))djAi{r) 

-edk, (A{k - A{r)) ■ F^,,{r, k - A{r)) + B{r) ■ M{k - A(r))) djAi{r) 

-sAi{k - A{r))(djdi(j>{r) - {WE{k - A{r)) x djB{r))^ 
+ sdjB{r)-M{k-A{r)), 

with the convention to sum over repeated indices. Substituting k = k — A{r) 
one obtains 

r-j = djE{k) -sd~^, (FLor(r, fc) • A{k) + B{r) ■ M(k)) 

and 

fcj. = kj-diAj{r)ri 

= -dj<f>{r)+diE{k)djAi{r) 

- e dk, (A{k) ■ Fj^orir, k) + M{k) ■ S(r)) djAi{r) 
+ eAi(k)dr^Fi^ori{r,k) + edjB{r) ■ M{k - A{r)) - diAj{r)fi 
= -d,^{r)+n[d,Aiir)-diA,{r)) 

+ eAi (k) dr^ FLor / (r, fc) + e a^- B (r ) • M (k) 
= - d,cl>{r) + (f X B{r))^ +eAi(k) dr^F^,n{r,k) +edjB{r) ■ M{k) , 
which, in more compact form, read 

r = VE{k)-eWi(A(k)-Fi^^r{r,k)+B{r)-M{k)^ , 

(59) 

k = -V(/)(r) + r X B{r) + eVr (A{k) ■ FLor(r, k) + B{r) ■ M(k)^ . 
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As the next step we perform the change of coordinates induced by T, 



q = r + eA{k) , p = k- A{r) + eV^ {A{k) ■ A{r)) , 



(60) 



and then switch to the kinetic momentum 



V 



P - Ail) 

k + eAi(k)\7Ai{r) ^ eAi(k.)diA{r) + 0{s^) 
k + eA(k) xB{r) + 0{e^), 



(61) 



where we used Taylor expansion. The inverse transformations are 



k = 



r 



q - e Aiv) + 0{e') , 
v-eA{v) X B{q) + 0{e^). 



Recall that we want to show that {q, v) satisfy the semiclassical equations of 
motion ij^l, where q is identified with r and v with k. The new notation is 
introduced here, only to make a clear distinction between the canonical variables 
(r, k) in the reference representation and the canonical variables (q, p) in the 
original representation. 

We now substitute (|6()|l and (|61|l . In the following computations we use 
several times Taylor expansion to first order and drop terms^ of order e^. In 
particular in the terms of order e one can replace r by g and fc by We find 



where it is used that v = i^Lor + 0(e). Thus we obtained the first equation of 
(O. For the second equation we find 



qj rj+eAj{v) 




- e5,^. (( - V0(q) + VE{v) X Biq))^Ai{v) + B{q) ■ M(«)) 



+ ediAjVi 



djE{v) - evi (djAi - diAj^ - e B{q) ■ djM{v) 
djE{v) - e{v X n{v)) . -eB{q) ■ djM{v) , 




- dj(f){q) + eAi{v)didjcf){q) 

+ {qx B{q))^ - e{A{v) X B{q))^ -e(qx {Ai{v)diB(q))) 



+ e Ai {v)dq^ FLori{q,v)+e djB{q) ■ M{v) 
+ e{A{v) X B{q))^+e(^A{v) X {qidiB{q))) 
- + {qx B{q))^ + e d,B{q) ■ M{v) , 
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where the term 

eA{v)(dq^Fi^orl{q,v)+didj4>{q)) =eAi{v)(q X d,B{q))^ + 0{e^) 

cancels the remainmg two terms. Changing back notation from {q,v) to (r, k), 
this concludes the proof of Theorem 1111 □ 

A Operator- valued Weyl calculus for r-equiva- 
riant symbols 

The pseudodifferential calculus for scalar-valued symbols defined on the phase 
space T*W^ = M^'^ can be translated to the phase space T*T<^ = T'^ x W^, T"^ 
a flat torus, by restricting to periodic functions and symbols. This approach is 
used by Gerard and Nier |GeNi| in the context of scattering theory in periodic 
media. 

In this appendix we present a similar approach to Weyl quantization of 
operator-valued symbols which are not exactly periodic, but T-equivariant with 
respect to some nontrivial representation r of the group of lattice translations. 
We obtain a pseudodifferential and semiclassical calculus which can be applied to 
r-equivariant symbols like the Schrodinger Hamiltonian with periodic potential 
in the Zak representation. In particular, the full computational power of the 
usual Weyl calculus is retained. The strategy is to use the strong results available 
for phase space K^'^ by restricting to functions which are r-equivariant in the 
configurational variable. 

Let r C M"* be a regular lattice generated through the basis {71, . . . ,7d}, 
7j e W^, i.e. 

r = e M'* : a; = Ej=i"j lo fo^' some a £ . 

Clearly the translations on by elements of F form an abelian group isomorphic 
to Z'*. The centered fundamental cell of F is denoted as 

M = {.T e M"^ : a; = Ejli^j Ij for a, e [-i, i]} . 

Let 7i be a separable Hilbert space and let r be a representation of F in 
B*{TL), the group of invertible elements of B{Ti.) , i.e. a group homomorphism 

T:T^B*{n), 7h^t(7). 

If more than one Hilbert space appears, then r denotes a collection of such 
representations, i.e. one on each Hilbert space. 

Warning: In the application of the results of this appendix to Bloch electrons 
the lattice F corresponds to the dual lattice F* in momentum space M.'^. 

Let be the operator of translation by 7 G F on 5(R'^, H), i.e. {L^ifi){x) = 
ip{x ~ 7), and extend it by duahty to distributions, i.e. for T e 5'(M'*,H) let 
(L^r)(^)=T(L_^^). 



28 



Definition 17. A tempered distribution T G S'{W^, H) is said to be T-equivari- 
ant if 

L^T = T{-y)T for all 7 G T , 

where {t{j)T){ip) = T(t(7)-V) for ip G S{R^,n). The subspace of r-equivari- 
ant distributions is denoted as <S^. Analogously we define 

= 1^ G Ll^{R'^, n) : il>{x - 7) = t(7) V(a;) for all 7 G t} , 

which, equipped with the inner product 

{'fi,ij)n^= dx{ifi{x),i){x))n , 
Jm 

is a Hilbert space. Clearly 

= {V' G C°°(M^H) : ^{x - 7) = r(7) for all 7 G r} , 

is a dense subspace of Hr- 

Notice that if r is a unitary representation, then for any G Hr the map 

X H-> {(p(x),ip{x))-H is periodic, since 

{(p{x - 7), ip{x - 7))-H = (t(7)(^(x), T(7)V'(a;))-H = (-^(x), ipix))n ■ 

Now that we have r-equivariant functions, we define r-equivariant symbols. To 
this end we first recall the definition of the standard symbol classes. 

Definition 18. A function w : R^'' [0, +00) is said to be an order function, 
if there exist constants Co > and A^o > such that 

w{x) <Co{x- y)^° w{y) 

for every x,y £ K^*^. 

It is obvious and will be used implicitly that the product of two order func- 
tions is again an order function. 

Definition 19. A function A G C°°{M.'^''',B{Hi,n2)) belongs to the symbol 
class S'^{B{Hi,H2)) with order function w, if for every a, /3 G N"^ there exists a 
positive constant Ca,0 such that 

\\{d^d^A){q,p)\\^^^^ ,^^^ < C„,p w{q,p) (62) 

for every g,p G R"^. 

Definition 20. A map A : [0, £0) Sw{B{'Hi,'H2)), £ ^ A^ is a scmiclassical 
symbol of order w, if there exists a sequence {Aj}j^^ c Aj G 5'™(S(Wi, W2)) 
such that 

00 

Ax^e-'A, in S^{B{ni,n2)), 
0=0 
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which means that for every n G N and for all a, /? G N"* there exists a constant 
Ca,i3,n such that for any e G [0, £o) one has 

ri-l 

d'^dP(A,{q,p)-Y,e^Aj{q,p)) <e"Cc,,f,,nw{q,p). (63) 

^ . _ / B('H^ .Tio) 



The space of semiclassical symbols of order w is denoted as S"^ {e, B{Hi,7i2)) 
or, if clear from the context or if no specification is required, as S^{e). The 
space of formal power series with coefficients in S'^{B{Ti,i,Ti,2)) is denoted as 
M^{e,B{Hi,H2)). 

Definition 21. A symbol G S^{e,B{'Hi,'H2)) is r-equivariant (more pre- 
cisely (ti, T2)-equivariant), if 

A{q -1,P) = 7-2(7) -4e((7,_p) ri(7)~^ for all 7 G F . 

The space of r-equivariant symbols is denoted as S'™(£, B(7ii, 7^2))- 

Notice that the coefficients in the asymptotic expansion of a r-equivariant 
semiclassical symbol must be as well r-equivariant, i.e. if A^ x X)j^o^"'^i' 
A, G S^{e,B{li.i,'H2)), then Aj G ^^(^^(Hi, H2)). 

Given any r-equivariant symbol A G 5^(i3(7ii, 7^2)), one can consider the 
usual Weyl quantization A, regarded as an operator acting on S'{W'-,TLi) with 
distributional integral kernel 

KA{x,y) = -^f d^A(i(:r + y),C) e'«(--^)/^ (64) 

Notice that integral kernel associated to a r-equivariant symbol A is r-equivari- 
ant in the following sense: 

Ka{x ~ '^,y - = T2{'y) KA{x,y)Ti{'jy^ for aU 7 G T . (65) 

The simple but important observation is that the space of r-equivariant 
distributions is invariant under the action of pseudodifferential operators with 
r-equivariant symbols. 

Proposition 22. Let A G Sl^{B{ni,n2)), then 

l5;^(R^7^i) C5;(M^H2). 

Proof. Since A maps S' {M.'^ ,Ti,i) continuously into iS'(M'^, 7^2), we only need 
to show that {L^AT){(p) = (r2(7)Ar)(v5) for aU T G 5;^(M'',Hi) and (p G 
S{R'^,n2). 

To this end notice that as acting on S{R'^,Ti,2) one finds by direct com- 
putation using (EH) that A* = (^(7)"^)* A* r2(7)*. Indeed, let -0 G 
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,'H2), then 

{A*Lj^){x) = I dyKA*ix,y)ij{y--/)^ I dy KA'{x,y + -f) ipiv) 

dy (n (7)"^)* {x-j,y) T2{j)* V'(y) 

Hence, using the fact that t is a representation and that L^T — Ti{"f)T, 
{L^AT){^) = T{A*L^^ip)=T{L^^nijrA*{T2{j)-'rip) 



□ 



iT2h)AT,{^)-'L^T){ip)^{T2h)AT){^). 



For the convenience of the reader we also recah the definition and the basic 
result about the Weyl product of semiclassical symbols. For a proof see e.g. 
El- 
Proposition 23. Let A e S""i(e,B(7^2,^3)) anrf S e S'™^ (e, H2)), then 
AB = C, with C G S'"^'"-'{e,B{Hi,H3)) given through 



C{e,q,p) = exp ( y (Vp ■ V, - V^. ■ V,) ) A{e,q,p)B{e,x,0 



=:AiB. 
(66) 



The corresponding product on the level of the formal power series is called 
Moyal product and denoted as 

tt : M^' (£, B(W2, W3)) X M"'^(e, BiHuH^)) ^ M^'^"'^ (e, B{nun^)) . 

The T-equivariance of symbols is preserved under the pointwise product, the 
Weyl product and the Moyal product. 

Proposition 24. Let A, e S^'^{e,B{n2/H:i)) and B, e 5^'^ (e, T^a)), 
then A,B, £ S^'"^^ {e, Bin^na)) andA.^B, e S'^'^'^'^ (e, T^s))- 

Proof. One has 

A{q-l,p)Be{q~l,p) = T3(7)Ae(g,p)T2(7)~V2(7)B£(g,p)Ti(7)-i 

T3{j)A,{q,p)B,{q,p)Ti{-f)-\ 

which shows Ai^B^ S S'^'i'^'^ (e, B(7ii, Tia)) and inserted into yields immedi- 
ately also A,Jb^ e S'^i'"2(e,i3(Hi,H3)). □ 

An analogous statement holds for the Moyal product of formal symbols. 
A not completely obvious fact is the following variant of the Calderon- 
Vaillancourt theorem. 
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Theorem 25. Let A G Sl{B{Ti.)) and ti,T2 unitary representations of T in 
Bin), then A £ B{Hr^,'Hr2) and for A^ £ Sl{e,B{n)) we have that 

sup ||1£||b(w^^,-h^^) < oo. 

66[0,6o) 

Proof. Fix n > d/2 and let w{x) = (x)^". We consider the weighted L^-space 



Let j — 1,2, then Ti^^ C and for any ip G Tirj one has the norm equivalence 

Ci UWn.^ < Uhi < C2 (67) 

for appropriate constants < Ci, C2 < 00. The first inequality in H67() is obvious 
and the second one follows by exploiting Tj-equivariance of "0 and unitarity of 



= J2 f dxw{xr\\T,i^r'ij{x)rn = j2 f dxw{xrmx)r^ 
< {^(^)'} / d^m^)\\H<c2\mH.^. 

According to H67|l it suffices to show that A G B{L^) and to estimate the norm 
of A^ in this space. 

Let ip G C^{W^,Ti.), then by the general theory Aip is smooth as well (see 

|Fo) . Corollary 2.62) and thus, according to Proposition 123 M G C^{R'^,H). 
Hence we can use H67II and find 

However, by Proposition l23l we have that w jj v4 (j w^^ G S^{e, B{Ti,)). Thus from 
the usual Calderon-Vaillancourt theorem it follows that 

Ik ^ '^~^lle(L2) - '-^d. II wM|w"^||^2d+l(^2d) • 

This shows that for A G S'i(^(H)) we have 1 G BiUr^Ur^)- Wiihwl A^lvj-^ G 
S^{e,B{n)) for A^ G Sl{e,B{rC)), we conclude that 

sup \\Ae\\B(H^^,H^^) < 00 
eG[0,£o) 

by the same argument. □ 

Remark 26. It is clear from the proof that the previous result still holds true 
under the weaker assumption that ri and T2 are uniformly bounded, i.e. that 

sup||Tj(7)lle(W) < C", j = l,2. 
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Finally we would also like to show that for A G S}^{B{TL)) the adjoint of 
A as an operator in B{Tlr): denoted by At, is given through the quantization 
of the pointwise adjoint, i.e. through A*. Here it is crucial that t is a unitary 
representation. 

Proposition 27. Let S^{B{Ti,)) with a unitary representation r (with ti = T2 = 
t) and let A'^ be the adjoint of A G B{TLr), then A'^ = A* . 

Proof. Let V e and 95 G such that ^ := Ia/ G C^(R'^, H), where 1a/ 
denotes the characteristic function of the set M . Such </? are dense in Tir and 
the corresponding ip can be used as a test function: 

= / dx{^{x),{Aij){x))^ dx{^{x),{A^){x)) 



dx{(A*^){x), ^{x))^ 

dx( dyKA-ix,y)lp{y),ip{x)) 



H 



dx( dyK\{x,y)ip{y),ij{x) 
^ Jm ' ^ 



/ XI ( / d2;i^l(.T + 7,2;)<y9(?;), ?/;(.T + 7) 
Jm \ Jm ' 

X ( / dyr"\7)i4:^(a;,?/ - 7)t(7)(^(?;), r"^(7)-!/;(a;)) 

/ ^^X\/ dy-ftri(a;,2/-7)(^(?/-7), -0(2;) 
Jm ^gp Wa/ ' 

dx/ / dyKA{x,y)ip{y),ip{x) 



M ' JR'' 

dx{(A^-p){x), tp{x))^ = {'^*^,i^)n 

M 

In particular, we used the r-equivariance of the kernel H65|l and of the functions 
in Tir and the unitarity of r. By density we have A* = A'' . □ 

B Hamiltonian formulation for the refined semi- 
classical model 

The dynamical equations ((Sj, which define the e-corrected semiclassical model, 
can be written as 

r = \7 ^Hscir, k) - e k X Qn{K.) , 

(68) 

k = —VrHsc{r, k) + r X B{r) 
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with 

H,,{r, k) := E^{k) + (/.(r) - eM„(/c) • B{r) . 

Recall that we are using the notation introduced in Remark [3 and that B and 
r2„ are the 2-forms corresponding to the magnetic field and to the curvature of 
the Berry connection, i.e. in components 

B{r),, = {d,A,-djA,){r) 

for i,j S {1, . . . , c?}, and 

VLn{n)^j ^ [d^Aj - djAi){K) . 

We fix the system of coordinates z = (r, k) in M.'^'^. The standard symplectic 
form 00 = Qoiz)im A dzi, where l,m € {!,..., 2d}, has coefficients given 
by the constant matrix 

where I is the identity matrix in Mat((i, R). The symplectic form, which turns 
H()8|) into Hamilton's equation of motion for Hsc, is given by the 2-form 03,6 ^ 
&B,e{z)im d^m A dz; with Coefficients 

For £ = the 2-form 0_B,e coincides with the magnetic symplectic form Qb 
usually employed to describe in a gauge-invariant way the motion of a particle 
in a magnetic field f |MaRa| . Section 6.6). For e small enough, the matrix H69|l 
defines a symplectic form, i.e. a closed non-degenerate 2-form. Indeed, since 
detOs = 1 it follows that, for e small enough, 03,6 is not degenerate. In 
particular it is sufficient to choose 

£< sup (||B(r) !]„(«:) II + ||^2„(^) II). 

The closedness of Qb,6 follows from the fact that B and fi^ correspond to closed 
2-forms over R'*. 

With these definitions the corresponding Hamiltonian equations are 

Qb,s{z) z = dHsc{z) , 

or equivalently 

{ B{r) -I \(r\^( VrH{r, k) \ 
\ I e ri„(«:) )\k ) \ V^H{r,K) ) ' 

which agrees with (|68l) . We notice that this discussion remains valid if ri„ admits 
a potential only locally, as it happens generically for magnetic Bloch bands. 
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